Evidence of fractional transport in point- vortex flow 



Xavier LeoncinO 

Courant Institute of Mathematical Sciences, New York University, 251 Mercer St., New York, NY 10012, USA 

Leonid Kuznetso\0| 

Lefschetz Center for Dynamical Systems, Division of Applied Mathematics, Brown University, Providence, RI 02912, USA 

George M. ZaslavskjOl 

Courant Institute of Mathematical Sciences, New York University, 
251 Mercer St., New York, NY 10012, USA and 
Department of Physics, New York University, 2-4 Washington Place, New York, NY 10003, USA 

Advection properties of passive particles in flows generated by point vortices are considered. 
Transport properties are anomalous with characteristic transport exponent fi ~ 1.5. This behavior is 
linked back to the presence of coherent fractal structures within the flow. A fractional kinetic analysis 
allows to link the characteristic transport exponent fi to the trapping time exponent 7 = 1 + ^. The 
quantitative agreement is found for different systems of vortices investigated and a clear signature 
is obtained of the fractional nature of transport in these flows. 



I. INTRODUCTION 

One area extensively studied in the last twenty years 
is the phenomenon of chaotic advection This phe- 

nomenon results from the chaotic nature of Lagrangian 
trajectories and enhances the mixing of tracers in lami- 
nar flows. Indeed in its absence the mixing relies on the 
less efficient molecular diffusion. Hence, its applications 
are important in geophysical flows where advected quan- 
tities vary from the ozone in the stratosphere to various 
pollutants in the atmosphere and ocean, or such scalar 
quantities as temperature or salinity. This growing in- 
terest in geophysical flows increases the relevance of two- 
dimensional models and more specifically the problem 
of advection in a system with many vortices |inl|-[T^. 
Moreover, different observations and numerous models 
have shown that the transport of advected particles is 
anomalous and can be linked to the Levy-type processes 
and their generalizations ^3|"|23- Another peculiarity 
of two-dimensional worlds is the presence of the inverse 
energy cascade in turbulent flows, which results in the 
emergence of coherent vortices, dominating the flow dy- 
namics IH-Iil- 

In order to tackle this problem and especially the 
anomalous features, our approach has been gradual and 
the present work completes a series of papers [3|-^3|, 
which consists of gradual successive steps of dynami- 
cal investigations of transport in two-dimensional flows. 
We restricted ourselves to a relatively simple model, and 
performed a thorough analysis of the dynamics of trac- 
ers. This approach allows to infer some properties of 
the kinetics which actually govern transport. To settle 
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for a model, we recall that systems point vortices have 
been used to mimic the dynamics of flnite-sized vortices 
[sBI-fsll, and mention that the evolution of 2D turbu- 
lence after the emergence of the vortices has been suc- 
cessfull y d escribed by punctuated Hamiltonian models 
[iMHSll^- Moreover, noisy point vortex dynamics have 
been recently used to describe exact unstationnary two- 
dimensional solution of the Navier-Stockes equation 
Therefore, despite their relative simplicity, it was clear 
that systems of point vortices managed to capture some 
of the essential features of two-dimensional flows, it there- 
fore seemed natural to consider these systems as the basis 
for our investigations. 

In the following we investigate the dynamics as well as 
the advection properties of point vortex systems, more 
speciflcally we address the problem with 3 point vortices 
as well as systems of 4 and 16 vortices. A system of three 
point-vortices is integrable and often generate periodic 
flows (in co-rotating reference frame! (4ll|-[4flj. We then 
can investigate the phase space of passive tracers with 
Poincare maps. A well-defined stochastic sea filled with 
islands of regular motion is observed, among these are 
special islands known as "vortex cores" surrounding each 
of the vortices. The non-uniformity of the phase space 
and the presence of islands of regular motion within the 
stochastic sea has a considerable impact on the transport 
properties of such systems. The phenomenon of stickiness 
on the boundaries of the islands generates strong "mem- 
ory effects" and transport is found to be anomalous. On 
the other hand, the motion of N point vortices on the 
plane is generically chaotic for > 4 46j-[48]. The peri- 
odicity is lost when considering a system of four vortices 
or more, but snapshots of the system have revealed the 
cores surrounding vortices are a robust feature (4fll|-[5(il|. 
the actual accessible phase space is in this sense non uni- 
form and stickiness around these cores has been observed 
[3^ . In fact, a refined analysis has revealed that in sys- 
tems of 4 vortices and more the cores are surrounded by 
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coherent jets, within which little dispersion occurs and 
thus which give rise to anomalous transport properties 

The goal of this paper is to propose a unified view 
of transport properties in these system of vortices. In 
Sec. m we start to recall briefly the general equations 
governing dynamics of point vortices as well as passive 
tracers. We then describe with more details the motion 
of vortices and tracers for some specific cases. We start 
with integrable systems of vortices namely, a system with 
three identical vortices, and a system of with two identi- 
cal vortices and one of the opposite sign, a system which 
allows to set up parameters in order to set the motion 
of the vortices on a collapse course. We then consider a 
system of four and sixteen identical vortices. In Sec. IIIII 
we investigate the transport properties of advected par- 
ticles in these systems and measure the typical transport 
exponent, finally in Sec. lIVI we develop on the fractional 
aspects of transport and Hnk its anomalous features to 
the fractal nature of the topology of the fiow. 



II. POINT- VORTEX AND PASSIVE TRACER 
MOTION 

A. Basic Equations 

Systems of point vortices are exact solutions of the 
two-dimensional Euler equation 

[f^,*] = (1) 



dt 



(2) 

where is the vorticity and \1/ is the stream function. 
The vortices describe the dynamics of the singular distri- 
bution of vorticity 
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where z locates a position in the complex plane, Za = 
Xa+Wa is the complex coordinate of the vortex a, and ka 
is its strength, in an ideal incompressible two-dimensional 
fiuid. This system can be described by a Hamiltonian of 
N interacting particles (see for instance referred to 
as a system of N point vortices. The system's evolution 
writes 

dH ^ „. dH 

(4) 
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where the couple {kaZa,Za) are the conjugate variables 
of the Hamiltonian H. The nature of the interaction de- 
pends on the geometry of the domain occupied by fiuid. 
For the case of an unbounded plane, the resulting com- 
plex velocity field v{z,t) at position z and time t is given 
by the sum of the individual vortex contributions: 



v{z,t) 



1 ^ 
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and the Hamiltonian becomes 

iJ = kakfj \n\za - 



(6) 
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The translational and rotational invariance of the Hamil- 
tonian H provides for the motion equations Q three 
other conserved quantities, besides the energy. 



TV 
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Among the different integrals of motion, there are three 
independent first integrals in involution: H , + 
and L^; consequently the motion of three vortices on the 
infinite plane is always integrable and chaos arises when 
N>A [43]. 

The evolution of a passive tracer is given by the advec- 
tion equation 



i = t) 



(8) 



where z{t) represent the position of the tracer at time i, 
and ti(z, t) is the velocity field Q. For a point vortex sys- 
tem, the velocity field is given by Eq. J^, and equation 
^ can be rewritten in a Hamiltonian form: 
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dz ' 
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dz 



where the stream function 
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^'(z,Z,t) = -— ^ ka\n\z - Za{t)\ 



(9) 



(10) 



acts as a Hamiltonian. The stream function depends on 
time through the vortex coordinates Za(t), implying a 
non-autonomous system. 

Due to chaotic nature of the evolutions we rely heav- 
ily on numerical simulations. The trajectories of the 
vortices as well as those of the passive tracers are in- 
tegrated numerically using the symplectic scheme de- 
scribed in .521 and which has already been successfully 
used in [laM El 111 IS 13 • 



B. 3- Vortex systems 

A general classification of different types of three vor- 
tex motion, as well as studies of special cases were ad- 
dressed by many authors [41. 43,. .44 .4,'Tj | . Among the dif- 
ferent possibilities quasiperiodic motion of the vortices is 
found generically for solutions for which the motion of the 
vortices is bound within a finite domain. We consider two 
different type systems. On one hand we consider a system 
with three identical vortices (see ^^), in this setting the 
two extreme regimes for the advection pattern of strong 
and weak chaos can be investigated with good control 
and analytical expression of vortex core is given. On the 



3 



other hand in order to stress our results we consider also 
3-vortex systems in the vicinity of a configuration leading 
to finite time singular solution. Indeed three vortices can 
be brought to a single point in finite time by their mu- 
tual interaction. Aref j42| points out that this result was 
already known to Grobli more than a century ago. This 
phenomenon, known as point-vortex collapse was studied 
in 41, 45, 53, 54). Thus, under certain conditions, which 
depends both on the initial conditions, and the vortex 
strengths, the motion is self similar, leading either to the 
collapse of the three vortices in a finite time, or by time 
reversal, to an infinite expansion of the triangle formed 
by the vortices. Let us re-derive these conditions of a 
collapse or an infinite expansion. For a system of three 
vortices in an unbounded domain, the invariance of the 
Hamiltonian © under translations allows a free choice 
of the coordinate origin, which we put to the center of 
vorticity (when it exists), the angular momentum in a 
frame independent form is then rewritten as: 



K 



- ^ kikj 



(11) 

(12) 



The first condition is immediate as for a collapse to occur 
the frame independent angular momentum ljll|l has to 
vanish. For the second condition we look for conditions 
resulting in a scale invariant Hamiltonian. With such 
spirit let us divide all lengths by a common factor A in 
the Hamiltonian H. We readily obtain H'{X) — H + 
{J2 kikj) In A, we then obtain the collapse conditions: 



K = 0, 



0, 



(13) 



j.e, the total angular momentum in its frame free form 
and the harmonic mean of the vortex strengths are both 
zero Ijl3|l . Near a collapse configuration, the two condi- 
tions lfT3jl allow two different ways to approach the sin- 
gularity. Namely, we can change initial conditions which 
changes the value of K, or change the vortex strength 
and modify the harmonic mean. These lead to different 
types of motion which have been classified and studied 
in 0. 

In all considered cases, the relative vortex motion is 
periodic, i.e. the vortex triangle repeats its shape after a 
time T. This does not imply a periodicity of the absolute 
motion, since the triangle is rotated by some angle 6 
during this time. In general, 8 is incommensurate with 
2tt, rendering a quasiperiodic time dependence of zi{t). 
Let us consider a reference frame, rotating around the 
center of vorticity with an angular velocity ft = Q/T. In 
this co-rotating reference frame, vortices return to their 
original positions in one period of relative motion T, their 
new coordinates 



are periodic functions of time. In the co-rotating frame 
the advection equation retains its Hamiltonian form with 
a new stream function VP. In this frame ^ is time-periodic 
(^'(z, z, t+T) — '^{z, z, t)) and well-developed techniques 
for periodically forced Hamiltonian systems can be used 
to study its solutions [l^. Note, that the one-period 
rotation angle 6 is defined modulo 2tt, making the choice 
of the co-rotating frame non-unique. 

Anomalous properties of tracer advection in a fiow gen- 
erated by the motion of three identical vortices were an- 
alyzed in j^^l- The vortices were initialized with such 
initial conditions creating a large region of chaotic tracer 
motion. The structure of the chaotic region is quite com- 
plex, with an infinite number of KAM-islands (stratified 
with regular trajectories) of different shapes and sizes 
embedded into it. To visualize these structures, we con- 
struct Poincare sections of tracer trajectories (in the co- 
rotating frame) . A Poincare section is defined as an orbit 
of a period-one (Poincare) map P 



Pzq = z{T, Zq) 
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z{T, Zq) 



(14) 
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where z{t, zq) denotes a solution z{t) with an initial con- 
dition z{t — 0) = Zq. Plots of Poincare sections for three 
identical vortices in the strong chaotic regime as well as 
for a system near a collapse configuration are shown in 
Figure . We recall that vortex and tracer trajecto- 
ries were computed using a symplectic Gauss-Legendre 
scheme The exact conservation of Poincare invari- 
ants by the symplectic scheme suppresses numerical dif- 
fusion, yielding high-resolution phase space portraits. 

The Poincare sections are presented in the Figure 
show an intricate mixture of regions with chaotic and 
regular tracer dynamics, typical for periodically forced 
Hamiltonian systems. All three phase portraits share 
common features with the advection patterns: the 
stochastic sea is bounded by a more or less circular do- 
main, there are a islands inside it, where the tracer's 
motion is predominantly regular. In particular, all three 
vortices are surrounded by robust near-circular islands, 
known as vortex cores. An expression of the radius of the 
cores is computed when the vortices are equal ^29], and 
the minimum inter- vortex distance through time provides 
as well a good upper estimate of the core radii (see for 
instance j^sli^l); this estimate proves also to be quanti- 
tative in a four vortex system [231 as well as in 16- vortex 
svstems^34j|. 



C. 4 and 16-vortex systems 

Due to the generic chaotic nature of 4-point vortex 
system, understanding the vortex motion necessitates a 
different approach than for integrable 3-vortex systems. 
For identical vortices it is possible to perform a canonical 
transformation of the vortex coordinates [i^l- For a 4- 
vortex system, this transformation results in an effective 
system with 2 degrees of freedom, providing a conceptu- 
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Figure 1: Poincare Map for strongly chaotic system (upper 
plot) and the system in the close to collapse configuration 
(lower plot). 



ally easier framework in which perform well defined two- 
dimensional Poincare sections can be performed 
To summarize the results obtained in |13,E3|, the motion 
is in general chaotic, except for some special initial condi- 
tions, for instance when the vortices are forming a square 
the motion is periodic and the vortices rotate on a circle, 
then symmetric deformation {z^ = —z\ and — —z-i) of 
the square lead to quasiperiodic motion (periodic motion 
in a given rotating frame) . 

As a prerequisite to our investigations on the advection 
of passive tracers, a basic understanding of the vortex 
subsystem behavior is necessary. For this matter, an ar- 
bitrary initial condition of the 4- vortex system is chosen. 



although the Poincare section is computed and the de- 
sired generic chaotic behavior verified jl^. As we evolve 
from the 4 vortex system to 16 vortices the phase space 
dimension is considerably increased and due to the long 
range interaction between vortices (see the Hamiltonian 
El the energy does not behave as an extensive variable. 
Thus, in order to keep some coherence between the four 
vortex system and the sixteen vortex one, we chose to 
keep the average area occupied by each vortex approx- 
imatively constant. The switch from 4 to 16 vortices 
can then be thought of as increasing the domain with 
non-zero vorticity while keeping the vorticity "uniform" 
within the patch. The initial condition is chosen ran- 
domly within a disk and there is no vortices with close 
neighborhood to avoid any possible forced pairing. After 
that all position are rescaled to match the condition of 
uniform vorticity. 

We shall now consider some specific behavior of these 
chaotic dynamics, namely the phenomenon of vortex 
pairing. The simulations performed in j^^, Ull indicate 
that long time vortex pairing exists, in fact the forma- 
tion of long-lived triplet (a system of 3 bound vortices) is 
also observed and thus the relevance of three vortex 
systems is confirmed. In fact, the formation of triplets 
or pair of vortices concentrates vorticity in small regions 
of the plane and in some sense is reminiscent of what is 
observed in 2D turbulence. However, since no quadru- 
ples (or even larger clusters) are observed, we may spec- 
ulate that this concentration of vorticity is intimately 
linked to the dynamics and long memory effects asso- 
ciated with stickiness. Indeed this phenomenon can be 
thought as sticking phenomenon to an object of lesser 
dimension than the actual phase space with some con- 
straints: the object is reached by generating subsystems 
(2 and 3 vortex systems) whose integrability is a good 
approximate for a fairly long time. For comparison we 
mention that a typical time of an eddy turnover used in 
[l^ . corresponds to a time of order At ~ 1 — 5 in these 
systems while triplets and pairing are typical times are 
at least 2 orders of magnitude larger. 

We have now a sufficient knowledge of the dynamics of 
the chaotic vortex systems, and move on to the behavior 
of passive tracers generated by these fiows. For the sys- 
tem of 4 point vortices, successive snapshots have shown 
that passive tracers can stick on the boundaries of cores 
and jump from one core to another core during a pairing 
or escape from the core due to perturbations . The fact 
that a tracer is able to escape from a core means that 
the surrounding regions of the cores are connected to the 
region of strong chaos. The results presented in [^^l in- 
dicate that these regions mix poorly with the region of 
strong chaos. One way to track this phenomenon is to use 
Finite-Time Lyapunov exponents (FTLE) and to elimi- 
nate domains of small values of the FTLE 
Once these exponents are measured from tracers' trajec- 
tories whose initial conditions are covering the plane, a 
scalar field distributed within the space of initial con- 
ditions is obtained and the two dimensional plot of the 
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Figure 2: Local snapshot of the system of 16 vortices with 
9. 10* tracers. The vortices are located with the "+" sign. We 
can see the cores surrounding the vortices. 
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Figure 3: Four consecutive snapshots for the four vortex sys- 
tem and 1000 particles. These correspond to four consecu- 
tive pairing. Tracers are initially placed around one vortex. 
As pairings occur, some jump from a vortex core to another 
where they remain after. While the vortex-pairing occurs 
some particles also escape from the cores. We notice that af- 
ter four pairings all cores have been "contaminated" and are 
populated with tracers originating from the first core, while 
about 10% of tracers have escaped from the sticky region sur- 
rounding the cores. 



scalar field reveals regions of vanishing FTLE, namely 
the cores surrounding the vortices and the far field re- 
gion. The cores are thus regions of small FTLE, meaning 
that two nearby trajectories are bound together for long 
times, and this despite the core's chaotic motion. These 
properties reveal typically a sharp change of the tracers 
dynamics as it crosses from the region of strong chaos to 



the core. 



III. TRANSPORT PROPERTIES 

A. Definitions 

In general the deterministic description of the motion 
of a passive particle in the chaotic region is impossible, a 
local instability produces exponential divergence of tra- 
jectories. Even the outcome of an idealized numerical 
experiment is non-deterministic, indeed in this situation 
round-off errors are creeping slowly but steadily from the 
smallest to the observable scale. The long-time behav- 
ior of tracer trajectories in the mixing region is therefore 
studied using a probabilistic approach. In the absence 
of long-term correlations, a kinetic description, which 
uses the Fokker-Plank-Kolmogorov equation and leads 
to Gaussian statistics, works fairly well for many 
situations. Yet in the present case, the topology of the 
advection pattern (Fig. and the trapping of tracers 
in the neighborhood of cores (see Fig. indicate that 
anomalous statistical properties of the tracers should be 
expected. The singular zones around KAM islands and 
the cores give often rise to a stickiness phenomenon and 
produce long-time correlations, which result in profound 
changes in the kinetics. In some cases these "memory ef- 
fects" can be accounted for by the modification of the dif- 
fusion coefficient in the FPK enua.tio n IfiMl^al. b ut often 
their influence is more profound 0, EaSH El , 
and leads to a super-diffusive behavior with faster than 
linear growth of the particle displacement variance: 



(15) 



the transport exponent exceeds the Gaussian value: 
/i > 1. 

For all considered case, the vortices are moving within 
a finite domain. It is thus important to define what quan- 
tities will be measured to characterize the transport prop- 
erties of the system. There has been evidence in |^ of 
radial diffusion, but the diffusion coefficient D is vanish- 
ing as i? ^ oo with typical behavior D ^ In the 
case of more than four vortices we can still expect a sim- 
ilar type of behavior. However the region far from the 
vortices is of little interest when one want to address the 
transport properties of typical geophysical fiows and the 
most relevant is the region being accessible to the vortices 
also called the region of "strong chaos". In the 3-vortex 
systems the question is even more crucial as the accessi- 
ble domain of tracers (the chaotic sea) reduces to a finite 
region surrounded by a KAM curve. One way around 
this finite domain problem is to focus our interest on the 
character of tracer rotation, and for that matter, we de- 
fine its azimuthal coordinate in the center of vorticity 
reference frame 



9{t) = Arg z 



(16) 
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to be a continuous function of time, i.e. 9{t) G {—oo, oo) 
keeps track of the number of revolutions performed by 
a tracer. However for many vortices we choose another 
quantity, namely we consider tracer transport by mea- 
suring the arclength s{t) of the path traveled by an indi- 
vidual tracer up to a time t. The arclength s{t) writes 



S^{t) = / V,{t')dt' , 

Jo 



(17) 



where Vi(t') is the absolute speed of the particle i at 
time t' . The main advantage of this quantity is that it 
is independent of the coordinate system and as such we 
can expect to infer intrinsic properties of the dynamics. 
The main observable characteristics will be moments of 
the angle 9{t) or distance s{t): 



Mg = {\x^{t) - {x^{t))\'' 



(18) 



where i corresponds either to the i-th vortex or a tracer 
in the field of 3, 4 or 16 vortices and x stands for s or 9. 
The averaging operator (• • • ) needs a special comment. 
Expecting anomalous transport one should be ready to 
have infinite moments starting from q > qo- To avoid any 
difficulty with infinite moments we consider truncated 
distribution function, which was discussed in details in 
[fill | and allows to satisfy the physical restriction of a finite 
velocity and and the finite time of our simulations. This 
condition actually put some constraints on the maximum 
meaningful value q* of q , and beyond q* the moments 
are basically monitoring the population of almost ballistic 
trajectories. 

Up to the mentioned constraints, we will always con- 
sider the operation of averaging to be performed over 
truncated distributions. In this perspective all moments 
are finite and one can expect 
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with, generally, fi{q) ^ q/2 as is expected from nor- 
mal diffusion. The nonlinear dependence of /i(g) is a 
signature of the multifractality of the transport, . For 
more information on the appearance of multi-fractal ki- 
netics and related transport see 62, 63]. Some authors 
use the notion of weak (/i(g) = const ■ q) and strong 
(m('Z) 7^ const ■ q) anomalous diffusion |^ or strong 
and weak self-similarity [fi3 |. 

Recently fHF, 'H^, two types of anomalous diffusion 
were distinguished by the behavior of the moments, 
and the notion of weak and strong anomalous diffusion 
[sEESl or identically strong and weak self-similarity 
introduced. When the evolution of all of the moments 
can be described by a single self-similarity exponent ly 
according to 



fi{q) =v -q 



(20) 



refers to "weak anomalous diffusion", whereas the case 
when v in ll2?Hl is not constant, i.e. 



is named "strong anomalous diffusion". This distinction 
is important since in the weak case the PDF must evolve 
in a self-similar way: 



P(x,t)=t--/(0, 



£,^t-^{x-{x)) 



(22) 



whereas non-constant v{q) in lf2T]l precludes such self- 
similarity (see also the discussions in for details 
about the non self-similar behavior). 

In all our simulations, the results show that the trans- 
port of tracers in point- vortex systems is strongly anoma- 
lous and super-diffusive, hence to avoid redundancy we 
graphically only present the transport properties of pas- 
sive tracers obtained for a 3-point vortex flow in Fig. ^ 
The behavior for the other considered systems is very 
similar, namely our results show, that for all considered 
systems ^i{q) is well approximated by a piecewise Hnear 
function of the form: 



vq 

q-c 



for 
if 



q<qc 
q> qc 



(23) 



where c is a constant, and q* is the crossover moment 
number. 

Even though we only considered few different initial 
condition for the different vortex system, it is reasonable 
to assume that the transport properties obtained for such 
systems are fairly general since all give the same kind 
anomalous behavior with a transport exponent more or 
less around /i(2) ± 1.5. 




Figure 4: The exponent ^(q) versus the moment order q for 

the angle distribution {{\e{t) -e{t)\'') ~ i^^'^) is plotted for a 

3-point vortex system {t > 1.5 10^). We notice two linear be- 

q<2 , ,i{q) = 0.75q strengths 
q>2 , fi{q) = IMq - Cte ' strengtns 

are (—0.3, 1, 1). The period of the motion is T = 17.53. 



haviors: 



/i(g) = qv{q) , iy{q) ^ const 



(21) 
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IV. FRACTIONAL ASPECTS OF TRANSPORT 

A. Poincare Recurrences and trapping times 
exponent 

The origin of the anomalous transport properties can 
be Hnked back to the intermittent character of tracer mo- 
tion. This phenomenon is characterized by an anomalous 
distribution of recurrences of the Poincare map of tracer 
trajectories for the three vortex systems, or the alge- 
braic decay of the density of trapping time within jets 
in many vortex systems. To define recurrences, we take 
a region B in the chaotic sea, and register all returns of 
a Poincare map trajectory into B. The length of a recur- 
rence is a time interval between two successive returns. 
The distribution of the recurrence times does not depend 
on the choice of the trajectory in the chaotic region, in 
dynamical systems with perfect mixing this distribution 
is Poissonian but the Hamiltonian systems with coexist- 
ing regular and chaotic motions exhibit a power-law tails 
in the distribution. Recurrence distributions for tracers 
in three vortex systems show that all distributions have 
long tails, indicating, that between the returns tracers are 
being trapped in long flights of highly correlated motion. 
Long recurrences are distributed according to a power 
law 

P(t) - T""^ (24) 

The measured values of the exponent 7 in Refs.Hilli 
are all around: 

7-2.5. (25) 

The value of the decay exponent 7 does not depend on 
the choice of the domain B as long as it is taken in the 
well-mixed region, away from the KAM-islands. 

In order to establish the origin of the long recurrences 
the corresponding Poincare cycles (i.e. parts of the map 
orbit between the returns) are plotted color-coded by 
the recurrence times r, see Fig. (only the cycles with 
T > 10'* » (r) « 850 are shown). The orbits of the long 
cycles are concentrated at the boundaries of the chaotic 
sea, around the islands of regular motion and at the 
external boundary. The longest recurrences (shown in 
red) correspond to the trajectories that penetrate deeply 
into the hierarchical island-around-island structures, see 
Fig. This phenomenon, known as stickiness (of KAM 
island boundaries) introduces long correlations to the 
tracer motion and leads to anomalous kinetic properties 
of chaotic trajectories. Another way to characterize this 
stickiness phenomenon using distribution of time aver- 
aged speed has been used in Js^, typically each sticky 
zone corresponds to a special time averaged speed dif- 
ferent than the whole phase space average, hence sticky 
portions of particles trajectories can be identified (see 
Fig.EJ. 

When dealing with more than three vortices the 
chaotic nature of the vortex systems does not allow to 
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Figure 5: Visualization of the fractal origin of anomalous 
transport in three vortex flows. The particles whose trajecto- 
ries have longest return times penetrate deeply into the hier- 
archical island-around-island structures (upper plot). Sticki- 
ness to islands corresponding to different type of characteris- 
tic average speed reveals the multi-fractal nature of transport 
(lower plot). 

use the Poincare map. However we can circumvent this 
problem by actually noticing that when dealing with 
more practical situations, we are facing a "coarse grained" 
phase space, and each point is actually a ball from which 
infinitely many real trajectories can depart. Given this 
fact we can imagine that two nearby real trajectories di- 
verge exponentially for a while but then get closer again 
without actually going to far from each other, a process 
which may happen over and over in the case of stickiness 
observed around islands in three vortex systems. From 
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Figure 6: Tail of the distribution of trapping time intervals 
. We notice a power-law decay, with some oscillations. 
Typical exponent is p{t) ~ t~'' with 7 ~ 2.823. 



ready illustrated in Fig. that the boundaries of the 
core exhibit the stickiness. However we insist on the fact 
that the stickiness of the system has been confirmed by 
looking for jets. In this sense the method used is rather 
general, and could be applied to other Hamiltonian sys- 
tems. 
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the "coarse grained" perspective those two real trajecto- 
ries are identical. We then can infer that there exists 
bunch of nearby trajectories which may remain within a 
given neighborhood for a given time, giving rise to what 
is commonly called a jet j^E^- Note that the stickiness 
to a randomly moving and not well determined in phase 
space coherent structure imposes existence of jets, while 
the opposite may not be the case. 

To actually measure the jets properties of the system, 
we use the following strategy. Let us consider a given 
trajectory r{t) evolving within the phase space. For each 
instant t, we consider a ball B{r{t), S) of radius S centered 
on our reference trajectory. We then start a number of 
trajectories within the ball at a given time, and measure 
the time it actually takes them to escape the ball, and 
have then access to trapping time distribution which is 
plotted for a system driven by four vortices in Fig. 
and log-log plot clearly shows the power-law decay of the 
trapping times, with typical exponent 7 = 2.82. 



B. Jets and stickiness 

Given the trapping time distribution it is possible to 
compute a distribution of finite-time Lyapunov exponent 
(see [34]), which identifies a threshold ao* and allows 
to dynamically detect a jet. Two different types of jets 
are identified, namely slow jets evolving in the region 
far from the vortices as well as fast jets localized on the 
boundaries of vortex cores. These last jets are profoundly 
influenced by the phenomenon of vortex pairing, and this 
phenomenon is at the origin of some mixing as it actually 
leads to the division and merging jets (see Fig. EJ. 

This localization of the jets conflrms the results al- 



Figure 7: Relative evolution of a tracer with respect to an- 
other within a long lived jets located in the far field region of 
the flow generated by four vortices. The distribution within 
the jets is gives rise to an hidden order organized as "ma- 
troshkas" (a nested set of jets with increasing radii). 

We now focus on the inner structure of the jet, a first 
plot of the evolution of a relative position of a tracer with 
respect to another tracer's position while both are within 
a jet. A plot of such structure is presented in Fig. [3 
where we discover a fine structure formed by a hierarchy 
of circular (tubular) jets within jets. This last structure 
found for a 4-vortex system, we also looked at a jet in a 
system of sixteen vortices. The results are plotted in Fig. 
IHl where the relative position of the ghost (relative tracer) 
is colored differently for different time periods of the life 
of the jet. We can see that effectively the nested set of jets 
within jets remains, and that the ghosts is also spiraHng 
back and forth in between. We also see the ghost going 
back very close to the tracer. This hierarchical structure 
is also reminiscent of the discrete renormalization group, 
and we speculate that log-periodic oscillation described 
in should be observed. 

C. Kinetics 

In some of the previous publications (see, for exam- 
ple, EM IZ^I) it was clearly indicated that the 
properties of anomalous transport are sensitive to phase 
space topology. More specifically, if we use the fractional 
kinetic equation [13; ES in the form 
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Figure 8: Jet structure for a long lived jet located in the 
region of strong chaos for a system of 16 vortices. The col- 
ors are characterizing different moments of the life of the jet 
corresponding to approximatively equal time intervals. They 
chronologically range as cyan<blue<green<red<magenta. 
We see a similar structure of jets within jets as observed in 
Fig. and the tracer spiraling back and forth between them. 



to describe distributions P{0,t) of rotations over angle 
9, then the transport coefficient V and exponents (a,/3) 
depend on the presence of different structures such as 
boundaries of the domain, islands, cantori, etc. 

The caption |3 shows that stickiness of trajectories to 
specific structures occurs with a filamentation of sticky 
domains along stable/unstable manifolds. In fact, dif- 
ferent sticky domains generate different intermittent sce- 
narios with some associated values of (a,/3) [6^ As 
a result, the real kinetics is multi-fractional and can be 
characterized by a set of values of (a, f3) or, more pre- 
cisely, by a spectral function of (a, /?) in the same sense as 
the spectral function for multi-fractals 0,123, Iz^l- The 
Figure characterizes the fact that trajectories, sticking 
to different structures (islands), have different angular 
velocities. Due to this, different asymptotic to the distri- 
bution function P{0, t) and different values of (a, /S) will 
appear for different time intervals. In other words, for 
a considered time interval one can expect a specific "in- 
termediate asymptotic" for P{0,t) and, correspondingly, 
different pairs {a, (3). Different classes of universality for 
the values {a,P) were discussed in We will now 

remind the consequences of some of these results. 

Multiplying Ij26|l by |6'|" and integrating it over 16*1 we 
obtain 

{\en ^ (27) 

or, in the case of self-similarity the transport exponent /i 



from the equation 



can be estimated as 



(28) 



(29) 



Expression II28II should be considered with some reserva- 
tions since the second and higher moment may diverge. 

It was shown in [23| that under special conditions the 
exponent 7 for the trapping time asymptotic distribu- 
tion 7/;(t) ^ T~'' can be linked to fractal time dimension. 
Moreover, 7 is related to the kinetic equation ll26ll as (fi^ 



/? = 7-l . 



(30) 



For the spatial distribution of particles, the simplest 
situation occurs when the diffusion process is Gaussian 
for which a — 2. When a hierarchical set of islands is 
present, a can be defined through scaling properties of 
the island areas. In the considered situation the random 
walk is more or less uniform in the regions where tra- 
jectories are entangled near stable/unstable manifolds. 
That means that we should expect that the value a ~ 2 
provides a good estimate. Finally, we arrive to: 



^i = 213 /a - 7 - 1 



(31) 



This last result is confirmed by all our observations for 
the different vortex systems studied in (scL l3^ 0| . We 
need to comment that it is not worthwhile to try to obtain 
/i with a high accuracy since a specific value of fx has no 
meaning due to multi-fractal nature of transport . 

In fact, we can be even more precise, as the exponent 
7 can be estimated to 7 « 2.5 leading to ^ « 1.5, which 
is a good approximation of the different observed values 
given the multifractality of the transport. We shall not 
reproduce here the estimation of 7 (see for instance (3^. 

for details) , but the idea revolves around the presence 
of an island of stability leading to baUistic or accelerator 
modes within the island. 



V. CONCLUSION 

In this paper the dynamical and statistical properties 
of the passive particle advection in fiows generated by 
three, four and sixteen point vortices has been reviewed. 
The goal of this work was to provide qualitative insights 
on general transport properties of two-dimensional fiows, 
more specifically geophysical ones, imposed by the topol- 
ogy of the phase space. The system of 16 vortices can be 
considered as a fairly large system while the 3- vortex one 
is the minimal one with non fixed distances between the 
vortices. Strong vortex-vortex correlation are observed. 
These correlations manifest themselves the formation of 
long lived pairs of vortices, triplets. Since these struc- 
tures are integrable, we can speculate that this form 
of stickiness occurs by forming quasi-integrable subsets. 
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The transport of passive tracers is found to be anoma- 
lous and super-diffusive in all investigated situation with 
a characteristic transport exponent ^ « 1.5. This phe- 
nomenon is explained by the existence of coherent jets , 
which are located in the sticky regions. 

These jets and sticky region exhibit a complex frac- 
tal structures in Figs El and |H1 and are responsible for 
the anomalous behavior. Indeed hence the power law be- 
havior observed in Fig. El characterize the exponent 7 
for trapping time, and the similarly power-law tails have 
been observed for Poincare recurrence times distribution 
in three vortex systems fsol Is^. In all situation we ob- 
serve an exceptional agreement with the 7 « /x(2) — 1 
relation resulting from the kinetic analysis performed in 
Sec. IIVI Hence since the notion of jet is quite general, we 



can say that the anomalous diffusion finds its origin in the 
existence of jets typically found around coherent struc- 
tures (islands and cores). Moreover we can access the 
transport properties of the global fiow by simply gather- 
ing escape time data from these coherent jets and using 
the equation (I^Tll resulting from the fractional kinetic 
equation lf2fijl . 

We therefore emphasize that the present work by ana- 
lyzing the role played in transport by the different fractal 
structures involved in our model-flows and by explicitly 
giving a typical value of the second moment exponent as 
well as the trapping time exponent should be of interest 
for the analysis of more realistic and complicated systems 
involving many vortices and coherent structures such as 
geophysical fluid dynamics. 
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